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The spatiotemporal structure formation problem is investigated in the region far above the trans- 
verse ferromagnetic resonance instability. The investigations are based on the dissipative Landau- 
Lifshitz equation and have been performed on a model which takes external fields, isotropic exchange 
fields, anisotropy fields and the demagnetizing part of the dipolar field into consideration. The nu- 
merical simulations for these models exhibit stationary domain structure in the rotating frame. 
Employing analytical methods and simplifying the model, certain features, such as the magnetiza- 
tion within the domains and the proportion of the system in each domain, are described analytically. 



PACS numbers: 76.50+g, 75.60.Ch 



I. INTRODUCTION 

Ferromagnetic systems driven by external magnetic 
fields have been under investigation for many years 0. 
Early work by Suhl already showed that there is a thresh- 
old of the pump amplitude where the uniform state be- 
comes unstable to homogeneous driving fields Q ■ When 
the difference between the pump field and its critical 
value is small (the weakly nonlinear case) , the experimen- 
tal findings can be analyzed by extensions of the theory 
of Suhl as reviewed in Q ■ The pattern formation in such 
a system can be described by amplitude equations 0, Q . 

However, when the probe is strongly driven, such per- 
turbative procedures break down. One-dimensional nu- 
merical simulations have predicted dynamic domains for 
a model including a transversely rotating field |fj. Dy- 
namic domains are stable solutions to the equations of 
motion. In the reference frame rotating with the driving 
field they exhibit a stationary domain structure, like that 
known from static domains, yet in the lab frame their 
position is stationary while the magnetization within the 
domains rotates at different angles. Plefka also obtained 
dynamic domains in one dimension by numerical simu- 
lation and was able to explain characteristic elements of 
the structure analytically Q for a very simplified model. 

The experimental difficulties involved with large pump 
amplitudes were solved already twenty years ago [8j • Re- 
newed interest in this topic has led to recent work, where, 
using the Faraday effect, dynamic domains were observed 
in garnet films driven by high power inhomogeneous driv- 
ing fields 

The length scales of the patterns we are interested in 
are large compared to the atomic distance of such sub- 
stances and so the relevant quantity is a macroscopic 
variable, the local magnetization m(r, t). Using a model 
that includes a saturating static field and a strong trans- 
verse pump field 0, we consider the structures occur- 



ring in a ferromagnetic film of side L with normal in the 
z-direction. 

The paper is organized as following: In Sec. II the 
physical model is described which enters in the Landau- 
Lifshitz Equation via the effective fields. The numerical 
results for the various dynamical domains are presented 
in Sec. Ill . Employing simplifications of the original 
model the linear stability of the structures found is con- 
sidered in Sec. IV leading to partial description of the 
numerical results. Conclusions are presented in Sec.V. 



II. LANDAU-LIFSHITZ EQUATION 

The dynamics of m(r, t) are described by the Landau- 
Lifshitz equation [lOllllj: 



dt m = m x h c ff — fm x (m x h e ff) 



(1) 
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r is a dimensionless damping coefficient. The effective 
magnetic field h c ff is made up of both external and in- 
ternal magnetic fields: 

h ff = He z + h (cos(ujt)e x + sin(u>t)e y ) (2) 
+ JV 2 m — m z e z + Km z e z . 

These terms are, in order of appearance: a static field 
of magnitude H in the z-direction (perpendicular to the 
plane of the film) ; an in-plane pump field with amplitude 
h and frequency u>; an isotropic exchange field (J > 0); 
the demagnetizing part of the dipolar field, describing 
the geometry of the film; and a uniaxial anisotropy, also 
in the z-direction. 

The explicit time dependence of is removed by a 
transformation to the rotating reference frame, yielding 

d t m = -m x (h cff - we 2 ) - Tm x (m x h eff ) , (3) 

while the effective magnetic field now has the form 

h e ff = He z + he x + JV 2 m — m z e z + Km z e z . (4) 
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FIG. 1: Dynamically stable domain states obtained at the 
same parameter values: w = H = 3, h = 2, J = 10, K = 5, 
F = 0.1. a) "spin-up" bubble, period L; b) horizontal stripes, 
period L/2; c) horizontal and vertical stripes, period L; d) 
diagonal stripes, period \J L/2; e) "spin-down" bubble, period 
L; f) "spin-down" bubble, period L/2 



An alternative, frequently used form of the Landau- 
Lifshitz equation is the Gilbert form 



dt m — I'm x dt m = 



(1 + r 2 )m x 



(5) 



i + r 2 



i + r 2 ' 



-m x e. 



This is obtained by multiplying equation @ from the 
left by Tm x and recalling that m • m = 1 . The damping 
coefficient T is the same in both the Gilbert form of the 
Landau-Lifshitz equation Jsj and the common form (J3J. 

Throughout this work we present results in the rotat- 
ing reference frame. We recall that static results in this 
frame of reference transform to structures rotating with 
frequency lo in the x-y plane in the laboratory frame. 
Hence a state that is statically stable in the rotating ref- 
erence frame is dynamically stable in the lab frame. 



III. NUMERICAL SIMULATION 

We simulate equation J3J) with Q on a vector-parallel 
computer and workstations, using a spectral approach 
to deal with the internal field contributions. The num- 
ber of mesh points for the two-dimensional system was 
varied up to 256 x 256 and it is seen that a minimum 
spatial grid size of 64 x 64 is necessary to recognize to 



good accuracy all magnetic structures that are also visi- 
ble with a larger number of mesh points. Periodic bound- 
ary conditions are used. For 64 x 64 mesh points, the 
boundary conditions also have no effect on the structures 
within the system. This is confirmed by also perform- 
ing the simulations with unpinned spins, and observing 
that the same results are obtained. A Euler integration 
scheme is sufficient and therefore used for the time inte- 
gration of the system. Indeed more advanced methods, 
like Runge-Kutta and NDSolve of Mathematica, give es- 
sentially identical numerical results. 

Throughout the course of a simulation we require that 
the magnitude of the magnetization remain constant 
m| = 1. However, the accumulation of numerical errors 
means that this is not the case, and the magnetization 
tends to drift away from the value one. Therefore we add 
an additional Bloch-likc damping term to the equation of 
motion JSJ that recalls the length of the magnetization 
vector back to one at each time step. This Bloch damp- 
ing, although a reaction to a numerical artifact, is also 
justified physically [Tlj . 

To investigate systematically the presence of dynam- 
ically stable structures, all numerical and physical pa- 
rameters except the amplitude of the driving field are 
kept constant (H = u = 3, J = 10, K = 5, T = 0.1), 
while h is varied from h — upwards. Depending on the 
initial states and on the value of h, different coexisting 
dynamically stable solutions arc found. 



A. Homogeneous State 

The solution found most frequently in numerical sim- 
ulations is the homogeneous or uniform solution. In this 
case, all the spins throughout the system rotate at the 
same frequency, in phase, and with the same constant 
m z . The dynamically stable homogeneous state is found 
numerically for values of h below a certain critical value 
h c i w 3.7 and above another critical value h C 2 ~ 4.9. 



B. Stripes and Bubbles 

Some of the two-domain dynamically stable structures 
also found are shown in Figure This is a shaded 
contour plot of the m z component of the magnetiza- 
tion throughout the system. The dark areas depicted are 
those where m z < while the light areas imply m z > 0. 
Figure Q makes no statement about the x and y com- 
ponents of the magnetization, which (apart from at the 
wall) are spatially homogeneous and vary periodically in 
time with period lo. This implies that the spins within 
the domains rotate in the x-y plane in phase and at the 
same frequency. 

The simulations were performed with a great many dif- 
ferent initial conditions. The results shown in Figure ^ 
are characteristic of all the dynamic domain results ob- 
tained, i.e. dynamically stable domain states either have 
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FIG. 2: Schematic view of where the homogeneous, stripe and 
bubble solutions are found numerically, h is the amplitude of 
the driving field 



either translational symmetry (b, c and d) or cylindri- 
cal symmetry (a, e and f) in the z-component of the 
magnetization. 

It is observed that dynamically stable domain states 
only exist between two values of the driving field, called 
h C 3 ss 0.2 and h c 4 «5. As we approach the upper bound 
h c 4 the domain states become ever more flattened out, 
until eventually at h c ± the domain state merges into the 
homogeneous state. This continuous transition is in di- 
rect contrast to the lower bound h C 3 where a well-defined 
domain state suddenly no longer exists. The values of 
h C 3 and h c i are seen to be the same whether bubbles or 
stripes are under investigation. 



C. Wall Structure 

The spatially distinct uniform regions are separated by 
a domain wall. Here the magnetization does not jump 
from the value in one domain to the next, but rather 
varies continuously across the finite wall width. The spins 
inside the wall rotate at the same frequency with the rest 
of the system, but with a phase shift compared to the 
spins within the domains. In the domain wall all the spins 
are parallel. In particular, this is also true for the bubble 
domain wall and so the whole structure including wall 
is not rotationally symmetric. The rotational symmetry 
of the system is broken by the presence of the driving 
field h. Such bubble domains are always obtained, even 
if rotationally symmetric initial conditions are used. 

The structure of the wall differs depending on h. For 
small h there is a large negative m y component of the 
magnetization. As h increases, m x dominates ever more 
at the center of the wall. In addition, the width of the 
wall increases with h. We noted above that the domains 
themselves become flattened out. This, together with 
the ever widening wall, means that it is difficult to say 
precisely when the domain state merges into the homo- 
geneous state at the upper bound h C 4- 



D. Hysteresis 

As indicated in Figure |2l there is tristability in the 
system, and which pattern actually occurs depends on 
the initial conditions of the system. Two hysteresis loops 
occurring in the system are shown in Figure [3] In partic- 
ular, the behavior of the system at h C 3 is of interest. We 



select as an initial state a domain state which would be 
stable just above h C 3. Reducing h a little below the crit- 
ical value h C 3, the two walls on cither side of the "spin- 
down" domain begin to move towards each other, and 
the magnetization within the walls to rotate faster and 
faster than the angular frequency of the driving field. 
Eventually the walls collide and leave the system in a 
homogeneous "spin-up" state. This process is shown in 
Figure The velocity of walls in this transient state is 
smaller the closer we are to the value h = h C 3, and it 
can take up to 60T, where T — 2n/u> is the period of the 
driving field, for the wall collision to occur. 



IV. ANALYTICAL INVESTIGATION 

A. Homogeneous Solution 

Homogeneous solutions must satisfy the following con- 
ditions 



V 2 m = 



m z 



and so the effective magnetic field becomes: 

h cff = He z + he x + (K - l)m z e z . (6) 

We insert © into equation and note that for steady 
state solutions {d^va = 0) the expression in square brack- 
ets in (J3J) must be parallel to m, with some unknown 
proportionality factor /x: 

he x + (S + (K - l)m z )e z + jm x e z = , (7) 

where we have introduced the following definitions: 



S = H- 



i + r 2 



7 



Tuj 

i + r 2 



(8) 



The parameter 6 is a measure of the deviation of the 
system from resonance and is called the detuning, while 
7 is the rescaled frequency. 

Taking the scalar product of equation Q with e x , e y 
and e z respectively and eliminating the unknown factor fi 
yields two equations for the 3 components of m. A third 
equation is obtained from the requirement that |m| = 1. 
We obtain a fourth order polynomial in m z : 



[(K - l) 2 + 7 2 ] m 4 z + 2S(K - l)m 3 z + [h 2 + 5 2 



21 4 

7 | m 

21 2 

7 rn z 



(K - If - 7 2 J m z z - 2S{K - l)m z - 5 Z = .(9) 

Either 2 or 4 solutions to equation (JHJ are found, depend- 
ing on the parameter regime. A saddle-node bifurcation 
separates the region of two solutions from that with four 
solutions. 

The stability of these four fixpoint solutions is dis- 
cussed by means of a linear stability analysis (see Ap- 
pendix). Figure [S] shows the components of the magneti- 
zation for that fixpoint where m z > 0, and the positions 
of h c \ and h C 2, between which this fixpoint is unstable to 
nonuniform perturbations with wavenumber k — > 0. For 
comparison, the dynamically stable homogeneous states 
found numerically are also shown. 
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FIG. 3: Hysteresis between the homogeneous and domain solutions. The figure plots m z . max — m z , m i n vs h, the amplitude of 
the driving field. For the homogeneous state m z .max TTLz, min = 0, while for bubble and stripe domains > 0. 

Each dot is the result of a numerical simulation. Arrows indicate the direction of change of the driving field h and the direction 
of transition. There is a sharp transition from the uniform state to domains at h c \ « 3.7 and from the domain state to the 
homogeneous state at /i C 3 ~ 0.2. The region around ft ps 5 is enlarged to show more clearly the sharp transition from the 
homogeneous state to the domain state at h C 2 » 4.9 and the continuous transition from the domain state to the homogeneous 
state at the somewhat higher value of /i C 4 
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t = 30T t = 40T 

FIG. 4: Disintegration of a spin-up bubble in time for a driv- 
ing field amplitude just below the critical amplitude h c z 

B. Domain Solution 

We make a simple ansatz to describe the magnetization 
in a two-domain state: 

ft , f mW < £ < Lq , . 

m( ^ = 1 mW < * < L(l - <?) (10) 

i.e. we consider only the magnetization within each do- 
main i = 1,2, and imagine ourselves to be a considerable 
distance away from the domain walls. The quantity q 
describes the proportion of the system in each domain. 

Again we look for steady state solutions (dtvay> = 0), 
and demand that the spatial inhomogeneity in each do- 
main vanish (V 2 mW = 0). The demand that the term 



FIG. 5: Dynamically stable homogeneous states at the pa- 
rameter values: H = lo — 3, J = 10, K = 5 for different 
values of h. The solid lines show the m x , m y and m z compo- 
nents of the magnetization determined analytically from Eq. 
©, while the dots represent results from numerical simula- 
tion. The critical values h c \ — 3.7 and h C 2 = 4.9 as shown are 
both determined analytically from Eq. iA7\ and observed in 
numerical simulations 

in square brackets in J5J be parallel to m yields: 

_ h W + ^_e 7 - -%m« xe 7 = ,i«m« (11) 

i + r 2 i + r 2 

with 

= He z + he x - m z e z + Km®e, . (12) 

These two equations for and m^ 2 ) are not indepen- 
dent, but rather are coupled via the demagnetizing term 
fn z which holds for the entire system and appears in each 
effective field h^l. Using the ansatz I|1U|I we write down 
the demagnetizing term as 

L 

m z = i j m z (0d£ = + (1 - q)mW . (13) 
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The other components of the magnetization are also sim- 
ply expressed as 



hK 



m. 



K 2 + 7 2 



in „ 



7/1 



K 2 + j 2 



(19) 



Figure shows the z-component of the magnetization 
within each domain, as given in l|18|) and compared to 
numerical simulations. We note that the domain solution 
can only exist for 



FIG. 6: The ^-component of the magnetization within dy- 
namic domains as obtained from numerical simulations and 
analytical calculation for different values of h. Each pair of 
dots for a given value of h corresponds to the m z values far 
from the domain wall for a single numerically stable simula- 
tion; the solid lines show the analytical result Eq. 1181 . The 
upper bound /i C 4 = 5.0037 is both calculated horn Eq. Il20ti 
and observed numerically, while the lower bound h c z = 0.2 is 
only observed numerically 



A further ansatz still needs to be made to obtain do- 
main solutions, and motivated by the symmetry of the 
numerical results we set: 



h < y/K 2 + 7 2 



(20) 



(1) 
m x 

■ 

m> y = 

(1) 

ml = 



(2) 
m x 

(2) 

m y 

(2) 



rn h 



(14) 



Although this symmetry is not inherent in the equations 
of motion, it is observed in all numerical simulations 
without exception and, more importantly, it fixes q and 
permits a simple closed solution form. 

By taking the scalar product of (|11|) for i = 1, 2 with 
e^, e y and e z and eliminating the unknown we find 
consistency only if 



= 5 



(15) 



This is an equilibrium condition for dynamic domains 
that fixes the proportion of the system in each domain. 
The effective magnetic field <|12|) in each domain i is now 



.0 

1 off 



(H - 5)e z + he x + Km^B z 



(16) 



Comparing equation (|16[) with JfjJ) , we see that all results 
from the homogeneous calculation in Section 11V Al can 
now be applied, providing H is replaced by H — S (thus 
S is replaced by (H - S) - = 0) and K by K + 1. 
Equation © then becomes 



to-' {K z +T 
with solutions 



'(h 2 -K 2 -7 2 )=0 (17) 



i«=±Wl- 



h 2 



K 2 



(18) 



A linear stability analysis (see Appendix) shows that the 
simple domain state l|10[l is stable for all values of h where 
the solution exists. This is in contrast to the numerical 
evidence, which indicates that there is a lower bound 
(shown as h = hc3 = 0.2 in Figure EJ) below which no 
numerical domain solutions exist. There is no indication 
of such a lower bound in any calculations carried out 
using the ansatz (|1C)|> . 

We recall that the ansatz (|10|) does not describe an en- 
tire dynamic domain, only that part far from the domain 
wall. Therefore any criteria for stability determined in 
performing a linear stability analysis on (|10fl will be only 
necessary conditions for stability, and not sufficient con- 
ditions. Critically, we have neglected any description of 
the shape or dynamics of the domain wall, and merely 
assumed a discontinuous transition from one domain to 
the next. In order to obtain a full description of dynamic 
domains, and, in particular, to determine the position of 
the lower bound h c ^, it will be necessary to include the 
wall and consider its effect on the stability of the system. 



V. CONCLUSIONS 

This work presents a model for a ferromagnetic film 
in which the exchange, dipolar and anisotropy effects are 
taken into account. Numerical simulation shows, depend- 
ing on the amplitude of the transverse pump field applied, 
different coexisting homogeneous states and dynamic do- 
main structures. The homogeneous state and its stability 
is wholly understood, and, using a simple model, certain 
features of the dynamic domain state are also clarified. 
The magnetization within the domains is determined, as 
is the proportion of the system in each dynamic domain. 
Nevertheless a lower bound for the existence of dynamic 
domains is observed in the numerical simulations that is 
not found with the simple domain-state model. 

An extension to this work has to include the effect of 
the domain wall, both its structure and its dynamics. 
Such an extension has been worked out [l2T | and will be 
published separately 01 leading to an understanding of 
the numerical results in their entirety. 
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APPENDIX A: LINEAR STABILITY ANALYSES 

1. Homogeneous Solution 

We denote the basic solution computed above as m . 
We then add a small time and space dependent pertur- 
bation to mo, in particular we consider one Fourier mode 
of such a perturbation: 



m(£, t) = m + <5m(£, t) = m + emi(t)e 



(Al) 



Now the magnetization m(£, t) must have magnitude 1. 
Thus mi must lie in a plane perpendicular to mo: 



m ( | • mi = 



(A2) 



We insert ansatz IjAlfl into the Landau-Lifshitz equation 
in the rotating reference frame @ and sort for powers of 
e. To order e: 

- rrii = mi x (h GfTi0 - we 2 ) + m x h cffj i + 

T [mi x (mo x h cff ,o) + m x (mi x h off , ) + 
m x (mo x h cff ,i)] , (A3) 

where we have introduced 

h e ff,o = He z + he x + {K - l)rn z . e z , 

h G ff,i = - Jk 2 mie lk ^ + Km z .ie z e lk( - - 8m z e z . 

The field h e ff,i depends on whether the case k = or 
k ^ is being considered, i.e. whether the perturbation 
<5m is uniform or non-uniform. If k = 0, 5m z = m z> \, 
whereas if k ^ 0, the fluctuations averaged over the entire 
system are zero and so Sm z — 0. Therefore 

for k = , h offj i = (K - l)m zA e z , 

for k ^ , h eff ,i = Km z ,ie z e iH - Jk 2 ime ik ^ . 

We introduce a coordinate system defined as 

ei = m , e 2 = m x e z , e 3 = m x (m x e 2 ) . (A4) 

In this coordinate system mi is written as 

mi(t) =/3(t)e 2 +7(t)e 3 

and equation (|A3J) assumes the form 

-/3e 2 - 7 e 3 = (h off , ■ m - wm Zi0 )(-/3e 3 + je 2 ) + 
r(m • h offi o)( / 9e 2 + 7e 3 ) + 
ei x [h cff ,i + Tei x h c ff,i] . (A5) 

We examine the four fixpoint solutions of © for stability 
with respect to uniform (k — 0) and nonuniform (k ^= 0) 



perturbations by considering the eigenvalues of the linear 
operator L defined by the right-hand side of equation 

The eigenvalues of the 2x2 matrix L for k = are 
positive for two of the four fixpoints. For the other two 
fixpoints, the trace of L for k ^ is negative for all values 
of fc, so we only need consider its determinant, given by 

Det(L kjt0 ) = (1 + T 2 ) (k 4 J 2 + k 2 JP + Q) , (A6) 
where 

P = (K(m 2 z - 1) + 2m • h offi0 ) - 2^m Z:0 , 
Q = (mo • h e ff,o) (m • h cff , + K(m 2 z - 1)) 
+pTO X)0 (m • h eff)0 + um a ,o) . 

Varying h, the determinant first becomes negative (and 
hence the fixpoint unstable) for k — ► at critical values 
of the driving field called h c \ and h C 2 given by 



P - \JP 2 - 4Q = 



Domain Solution 



(A7) 



The stability analysis of the domain solution is per- 
formed analogously to that of the homogeneous solution 
in Section IA II A small perturbation is added to the 
magnetization in each domain i = 1,2. As the Fourier 
modes in the first domain would couple to all of those in 
the second domain, we instead perform the calculation 
in real space, writing 



m«(e,i)=/3(e,<) W e^+ 7 (C,t)^e ; 

We obtain a system of 2 equations for each (/3^,7^), 
equivalent to (|A5(I 



(0 D (0 



-/3«e«- 7 «e« = 



m rl eff,0 



r(m«.h« !0 )(/3We«+ 7 (l) e«) (A8) 
fK |0 -l) 7 (i '-(m z )J (e«+re«) 



+J 



/3"W(e« _r e W)+ 7 " (i) (- 



-Te 



where h^g is defined in l|16|l . 

In the case of uniform perturbations, equation (|A8|) 
reduces to a system of four coupled ordinary differential 
equations. Numerically computing the eigenvalues of this 
4x4 matrix shows that all eigenvalues are negative for 
h < h C 4 = y/ K 2 + 7 2 and so the domain state is stable 
to (local) uniform perturbations. 

For non-uniform perturbations, we descritize equation 
(IA8|) by replacing the spatial derivatives with finite differ- 
ence expressions. For all values of h < h C 4 — \JK 2 + ■y 2 , 
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the real part of the eigenvalue spectrum is strictly neg- 
ative for both solutions mi in 1)18(1 . Hence the domain 
state is also stable to (local) non-uniform perturbations. 

Finally, we may consider what happens when the posi- 
tion of the wall is displaced from its equilibrium position 
given by (fl5|) . i.e., additionally letting q —> qo + eqi(t). 
Again, we obtain a linearized equation of motion. Inte- 
grating this equation over a small section of the contain- 
ing the wall, and then letting the size of this section go 
to zero 



Lqo+e 



lim 



d£ 



we obtain a single equation for qi(t): 



Qi = ~K [1 - T(l - m Zj0 )] qi 



(A9) 



The expression in square brackets in (|A9(I is strictly pos- 
itive, and so again the domain solution is stable to dis- 
placement of the wall from its equilibrium position. 
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